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O Abstract 

(N 

^ This book chapter gives an introductory review of a numerical framework for finding 

^ static and stationary vacuum black hole solutions. Such methods may be applied 

)—i to explore the exotic space of higher dimensional black holes that are thought to 

p exist. 

(N 

^ . 1 Introduction 

^ Whilst the black holes of 4D are well mannered, being spherically symmet- 

, ^^ ric or having special algebraic properties which enables them to be found 

analytically, moving beyond 4D many solutions of interest appear to have 
K^ no manners whatsoever. The problem of finding these unruly black holes be- 

(y^ comes that of solving the non-linear coupled set of partial differential equa- 

tions (PDEs) for the metric components given by the Einstein equations. 
In general it is unlikely that closed form analytic solutions will be found 
for many of the exotic black holes discussed earlier in this book. If we are 
(3 to understand their properties then we must turn to numerical techniques 

^~~' to tackle these PDEs. It is the purpose of this chapter to develop general 

• • numerical methods to address the problem of finding static and stationary 

. ;_H black holes. 

/\ Surely the phrase 'the devil is in the details' could not have a truer ap- 

^ plication than to numerics. The emphasis of this chapter will be to provide 

a road map to tackling the problem, where we formulate the problem in as 
unified, elegant and geometric a way as possible. We will also discuss con- 
crete algorithms to solve the resulting problem, but the extensive details of 
implementation will not be addressed, probably much to the readers relief. 
Such details can be found in the various articles cited in this chapter. 
Given the extensive efforts that have gone into understanding dynamical 

Chapter of the book Black Holes in Higher Dimensions to be published by Cambridge University 
Press (editor: G. Horowitz) 
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numerical simulation of gravity, the most obvious approach to find static 
and stationary solutions is to simulate a dynamical collapse of matter that 
is likely to form a solution of interest, or alternatively to simulate vacuum 
gravity starting with initial data that appears to already contain the black 
hole. Such a dynamical approach is indeed possible but is not the approach 
we will develop here for several reasons. Firstly, whilst such dynamical evolu- 
tions are well understood they are very complicated, typically requiring large 
computer resources. An important point is that in order to find an accurate 
stationary solution one would have to run a dynamical simulation for a long 
time to ensure that the resulting solution has indeed settled down, losing all 
its excitations to gravity waves, and this can be a serious challenge. Secondly, 
since one necessarily would be following the detailed evolution to form a hori- 
zon, to watch it radiate gravity waves and ring down, one is certainly doing 
far more work than is necessary if all that is required is the final stationary 
solution. Thirdly, many of the families of exotic solutions that have been 
discussed are likely to be unstable, at least for certain ranges of parameters. 
It is still interesting to find such black holes in order to elucidate the global 
structure of the space of solutions. For example in Chapter 4 we discussed 
that in Kaluza-Klein theory the inhomogeneous black strings and some of 
the localised black holes are thought to be unstable. However to confirm 
the elegant picture that these solutions are continuously connected to each 
other requires us to find these unstable solutions. In principle one might try 
to tune initial data to find an unstable solution, although in practice this is 
likely to be very hard. We also mention a more speculative concern, namely 
that in higher dimensions it appears that cosmic censorship does not hold 
[1] , and therefore it is unclear how generically one might expect to encounter 
singularity formation in dynamics. 

For all these reasons here we focus on finding static or stationary black 
holes by solving the static or stationary Einstein equations directly. Since 
these are PDEs the most important question is what character they have, 
as this will determine how to approach the problem. In the context of dy- 
namics the Einstein equations should be thought of as having hyperbolic 
character i.e. on small scales one has wave propagation along a light-cone. 
Consequently the dynamical Einstein equations are thought of as an ini- 
tial value problem, with data being specified on a Cauchy surface or past 
light cone. In contrast to this the static and stationary problem should be 
thought of as having elliptic character. One solves elliptic systems as bound- 
ary value problems where for a second order elliptic system one piece of data, 
for example Dirichlet/Neumann/Robin/oblique, is given on all boundaries. 
Physically in our static or stationary context such boundary conditions will 
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correspond to ensuring horizon regularity whilst also prescribing some par- 
ticular asymptotic behaviour. 

Our focus is to present the problem in an elliptic framework where one 
requires only entirely conventional techniques and modest desktop comput- 
ing resources. In particular we will discuss standard methods to solve the 
elliptic systems, namely relaxation and Newton's method. In addition the 
approach is amenable to the full variety of methods to represent the solu- 
tions, for example the classical finite difference approach, or more modern 
spectral, pseudo-spectral and finite element methods. 

Before embarking on our discussion of the numerical stationary problem 
for higher dimensional black holes some comments on the history of this 
problem are in order. Only with the rather recent revelation that black hole 
uniqueness breaks down in dimensions D > A did numerical exploration be- 
gin of gravity in higher dimensions. Of course in 4 dimensions uniqueness of 
Kerr meant that numerical work was not traditionally directed at vacuum 
black holes. Indeed a key ingredient of the classical proofs of uniqueness was 
to formulate the 4D stationary axisymmetric vacuum problem as an elliptic 
system [2]. However the numerical stationary elliptic problem has featured 
prominantly in the context of relativistic stars, particularly in the axisym- 
metric context, and there is a distinguished history of numerical work in 
this area beginning in the early 1970's (see for example the seminal works 
[21 m E] ) • For a review of this fascinating field the reader is referred to the 
Living Review article [5] . In the context of non- vacuum static and station- 
ary black holes such numerical methods were applied in [7j and p] to find 
exotic 4D charged solutions. In all these cases the approaches developed are 
generally based on the Weyl-Papapetrou form of the metric, exploiting the 
fact that the metric only depends non-trivially on two coordinates and has 
two commuting Killing vectors, so that solving the Einstein equations is 
consistently reduced to an elliptic problem. In the analytic context recent 
uniqueness theorems for vacuum stationary black holes in higher dimensions 
PlfTDj also phrase the problem as an elliptic system, and so far have restricted 
attention to D dimensional spacetimes with D-2 commuting Killing vectors, 
where again the Weyl-Papapetrou form is used (see also [llj). 

Numerical methods for metrics depending non-trivally on two coordinates, 
but without the restriction of having D-2 commuting Killing vectors were 
given in |12i [13] and applied to higher dimensional black holes. An analog 
of the Weyl form was employed for the metric. It was shown that a subset 
of the Einstein equations are elliptic in the metric components and may 
be solved regarding the remaining Einstein equations as constraints. One 
must show these constraints may be consistently solved by consideration 
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of the various boundary conditions. Being based on this analog Weyl form, 
these methods are manifestly non-covariant and by construction can only 
be applied to problems depending non-trivially on two coordinates. Whilst 
they have yielded interesting results they give a rather unstable numerical 
scheme for solutions in D > 4 which have axes of rotational symmetry 
(such as localised Kaluza-Klein black holes), presumably due to the lack of 
covariance, and are therefore hard to use in practice. We emphasize that the 
methods we will develop here will be fully covariant, and may be applied 
to problems depending non-trivially on an arbitrary number of coordinates. 
In our experience they are much better behaved in practice, for example, 
having no difficulties with axes of symmetry. 

For simplicity this chapter will focus entirely on finding vacuum black 
holes with no cosmological term and with a single component non-extremal 
Killing horizon. We will spend considerable time treating the case of static 
black holes, where most previous work has been directed. Indeed the numer- 
ical solutions of [11] discussed in Chapter 4 where found precisely using the 
methods we discuss. By analytic continuation of time this problem has an 
elegant geometric formulation as finding Ricci flat Riemannian solutions. We 
will study in detail the issue of formulation as an elliptic boundary value 
problem and give two algorithms to solve the resulting PDEs. In the re- 
mainder of the chapter we will show that the stationary case, which must 
directly be treated in Lorentzian signature, can also be thought of as an 
elliptic boundary value problem. 



2 Static vacuum black holes 

In this section we treat the static vacuum case following the approach of [U] 
with Headrick and Kitchen. Let us consider a general non-extremal static 
black hole solution with a single component horizon so that we may write 
the metric as, 

ds'^ = -N{xfdt^ + hij{x)dx'dx^ (2.1) 

where d/dt is the static timelike Killing vector, and at the horizon the norm 
of this vector vanishes so iV = 0. The zeroth law implies that the surface 
gravity given by the function k. = dnN\i\i=Q, where n is the unit normal 
vector to the horizon in a constant t slice, is actually a constant. A standard 
result of Euclidean quantum gravity is that any such static black hole may 
be analytically continued to imaginary time t = it to yield a Riemannian 
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manifold with metric, 

ds'^ = +N{xfdT^ + hij{x)dx'dx^ (2.2) 

where upon making r an angular coordinate with period t r-^ t + 2it/k. the 
metric at the horizon is smooth with no boundary there. To manifest this 
we may take Gaussian normal coordinates to the horizon where x* = {r, x"} 
and the horizon is located at r = 0. Then near the horizon the metric goes 
as, 

ds'^ ~ {K^^dT^ + dr^) + hab{r, x)dx''dx^ (2.3) 

and we see that the shrinking Euclidean time circle forms the angle of polar 
coordinates in M^, with r the radial coordinate. Whilst this polar coordinate 
system breaks down at the origin r = 0, simply by taking 'Cartesian' coor- 
dinates X = r cos KT and Y = r sin kt, one can write the metric in a chart 
that covers the horizon. 

Thus a static black hole can be written as a smooth Euclidean geometry, 
where Euclidean time r is periodic, and d/dr generates a C/(l) isometry. 
There is no boundary at the horizon and the geometry is perfectly smooth 
there. The horizon is marked by the vanishing of the vector d/dr, and con- 
sequently the horizon forms the fixed point set of the static isometry. 

Solving the vacuum Einstein equation without cosmological term is equiv- 
alent to finding a geometry that is Ricci flat, so R^i, = 0. Finding a vacuum 
static black hole solution can then be viewed as part of the more general 
problem of finding Ricci flat Riemannian geometries, the only special fea- 
ture of a static black hole geometry being the U{1) isometry, generated by 
the hyper-surface orthogonal vector d/dr, which leaves a codimension two 
submanifold (the horizon) fixed. 

There are several attractive features of this way of thinking. Firstly there 
is in principle no boundary at the horizon, and from the point of view of 
our boundary value problem, no boundary condition to impose there. The 
'in principle' qualification, which we will explain in more depth later, refers 
to the fact that if one wishes to use coordinates that manifest the U{1) 
isometry then these will not cover the horizon, in exactly the same way 
as polar coordinates do not cover the origin point. Whilst in practice it is 
sensible to adapt coordinates to any isometries and we shall deal with this 
in detail later, for the time being let us take a more formal view ignoring 
issues of implementation such as which coordinates to take. In principle one 
can always take coordinates that are regular at the horizon (the 'Cartesian' 
coordinates above) and then there is no boundary. 
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The second attractive feature is that the Euchdean continuation is pre- 
cisely what is done in semi-classical quantum gravity to consider the canoni- 
cal ensemble i.e. to work at fixed finite temperature. In particular the proper 
size of the Euclidean time circle asymptotically has the interpretation of 
being the inverse temperature of the black hole. Far from the black hole 
horizon we will impose boundary conditions on the metric that involve fix- 
ing amongst other things this proper size of the time circle. This boundary 
value formulation then naturally leads us to fix a very physical quantity, the 
temperature. 

Thirdly, the problem of finding Ricci flat Riemannian metrics is one with 
wider application than just finding static black holes. For example, numeri- 
cally finding the exotic Calabi-Yau geometries |15l [T6] that underpin certain 
string theory compactifications is currently required to make detailed pre- 
dictions of low energy phenomenology. In addition, particularly in the case 
of such Kahler metrics there may be new ways to think about the problem 
inspired by the geometry [n\ [T8] which might in the future be useful in our 
black hole context. 

Let us now consider how to formulate the problem of finding Ricci flat 
Riemannian geometries as an elliptic boundary value problem. 



2.1 The Harmonic Einstein equation 

The vacuum Einstein equation i^^j^ = is a second order quasi-linear PDE 
in the metric components. Perturb the metric about some background g^i, 
by a perturbation h^u, and then, 

where 

AlV = -2^^V-^A^^«a+%''/Iz.)k, Vf^ = yuh%--df,h, (2.5) 

and Al is the usual Lichnerowicz operator. The principle part of Ar, which 
we denote Pg, is given locally by taking only the two derivative terms, 

Pghf,u = - (g°''^df,dahi3ij + g"'^dudahi3^ - g'^'^dadph^u - g°''^d^duK^ 

(2.6) 
and this linear operator controls the very short wavelength behaviour of 
perturbations which determines the character of the equations R^^ about the 
background g. The condition that Rfj,^ = is elliptic about some background 
g gives the requirement that if one takes h^^y = a^j^e*'^"^'" for some constants 
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a^jy and any real non-zero A;^, then Pgh^y 7^ everywhere. Physicahy this 
condition means that nowhere can we find a point where short wavelength 
perturbations in a particular direction propagate as a wave. 

We see that for perturbations of the form h^y = 9(^m^) , where u is some 
vector field, then Pgh^v = 0. Thus the Ricci flatness condition R^u = is not 
an elliptic equation. Such a perturbation can be thought of as a short wave- 
length infinitesimal diffeomorphism generated by u. This diffeomorphism is 
given by h^y = Vr^Uj^) but for a u which varies on very short scales, then 
V(^tii/) ~ d(^^Ui^y Hence we may see lack of ellipticity of R^i, = as a conse- 
quence of gauge invariance. Without an elliptic set of PDEs we cannot treat 
the system as a boundary value problem as we wish to. In order to proceed 
we must therefore break this gauge invariance. However lifting gauge invari- 
ance need not imply breaking covariance. We emphasise that the method 
described below will indeed be fully covariant. 

Instead of considering the vacuum Einstein equation i?^,^ = 0, we will 
consider what we term the Harmonic Einstein equationPi?^ = where, 

< = R,.. - V(^^.) , r = g'-" (rv - f V) . (2.7) 

r is our usual Levi-Civita connection of 5, and T is another connection which 
we are free to choose, and then consider fixed. We term f the reference 
connection. Being constructed from the difference of two connections, the 
quantity ^ is a globally defined vector field. The equations Rj^^^ = have the 
great virtue of being elliptic. The principle part of the linearisation about a 
background g is simply, 

Pf V = -la'^^dc^dph^u (2.8) 

and thus for any Riemannian background g this clearly is elliptic. Taking 
hfj,i> = a^j^e*'^"^", then P^^h^y = a^yk'^ka which indeed only vanishes for 
vanishing /i^j^ or /c^ as required for ellipticity. An important point is that 
we have used the fact that we have analytically continued to Euclidean 
signature. For a Lorentzian signature metric g, then P^h^u = at a point 
X if we pick any non-zero but null vector k. In this Lorentzian case the 
Harmonic Einstein equation has hyperbolic character. 

The Harmonic Einstein equation has been used in the context of anal- 
ysis for decades for both the Riemannian elliptic problem as well as the 
Lorentzian dynamical hyperbolic problem (see for example [ID])) with varia- 

^ Note that we have referred to this equation also as the Einstein-DeTurck equation in the work 

[llllii. 
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tions in the precise definition of tlie vector that ah lead to elhptic/hyperbolic 
equations at least in the neighbourhood of a Ricci flat solution. 



The choice of vector field we have taken above in (2.7) is due to DeTurck 
who introduced it in the Riemannian context and later used it to show Ricci 
flow is parabolic as we discuss later [21]. For simplicity in what follows we 
reduce the freedom in the definition above and will take T to be the Levi- 
Civita connection of a reference metric g which we are free to choose and 
then consider fixed [22] . In this case we may write, 

Cm = 5"^ (V(,5^)^ - ^V^gap] , (2.9) 

where V is the covariant derivative of the metric g. q] 

In D dimensions there are D local coordinate degrees of freedom to fix in 
order to lift the gauge invariance of the Ricci fiatness condition. The con- 
dition ^'^ = precisely provides these D additional local conditions. The 
DeTurck choice of ^ can be thought of as a global version of the generalised 
harmonic coordinates introduced by Friedrich in the Lorentzian hyperbolic 
context and independently employed later by Garfinkle for numerical evo- 
lutions |24| [25] and used extensively since then. In generalised harmonic 
coordinates one writes ^° = fi^'^^F" + H"^ in local coordinates for some 
choice of H", although we note that H°' is not a vector field globally. Lo- 
cally this is the same as the definition above, where H" = —g^'^T". The 
vanishing of ^^ can be thought of as a generalised harmonic gauge condition. 
^ Given a chart the coordinates x° are functions over the part of the mani- 
fold covered by the chart. If in this chart we took F" = 0, then vanishing 
^ would imply V|x" = where V| is the scalar Laplacian, and hence the 
local coordinates would be harmonic functions - so called harmonic coor- 
dinates. For a general choice of reference connection, we have V^x" = H° 
corresponding to generalised harmonic coordinates. From now on we shall 



El 



refer to a metric with ^ = for our DeTurck choice of ^ in (2.7) as being in 
the generalised harmonic gauge. £j 

^ We note this is close to, although not the same as the Bianchi choice of vector field used, for 
exam ple, in 2Z . In particular the Bianchi choice does not lead to the simple principle symbol 
| |2.8| l except for metrics g close to the reference metric g. 

^ Note the subtle difference between the usual use of generalised harmonic coordinates, where 
one locally fixes H, and the DeTurck case here where instead one fixes f. These are 
inequivalent as the relation between them involves the metric g. 

* Since § = can be viewed as elliptic equations for the coordinate functions we see that whilst 
it is able to constrain the local degrees of freedom in the gauge, as with any elliptic equation 
one must specify boundary data, and this global gauge freedom must be fixed. When one 
specifies charts, one must specify the domain of these charts in M^ and this data precisely 
gives Dirichlet boundary conditions for the harmonic coordinate functions of the chart. This 
Dirichlet data for all the coordinate functions then implies all global data specifying the gauge 
is used up. Thus in the elliptic context harmonic coordinates are convenient as you have 
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A Ricci flat solution in generalised harmonic gauge ^ = does solve the 
Harmonic Einstein equation i?^ = 0. However the careful reader will notice 
that whilst the PDEs i?^ = are indeed elliptic for a Riemannian g (or 
indeed hyperbolic for Lorentzian g) there is no reason to suspect a priori 
that a solution to i?^ = has anything to do with a solution to the Ricci 
flatness condition. 

We now consider why solving i?^ = might lead to a Ricci flat solution 
presented in generalised harmonic coordinates. The situation is simplest in 
the Lorentzian hyperbolic context where the answer lies in the contracted 
Bianchi identity applied to the Harmonic Einstein equation. In either signa- 
ture this yields the linear PDE 

V'eM + V^- = (2.10) 

for the vector ^^. If one ensures ^^ and its normal derivative vanish on a 
Cauchy surface, then since the linear equation above is a wave equation, i.e. 
hyperbolic, for Lorentzian g then ^^ must remain zero under evolution of 
the metric in time. Hence in the dynamical hyperbolic context one simply 
imposes the vanishing of S, and its time derivative as constraints on the 
initial data for the metric, then solves the hyperbolic Harmonic Einstein 
equation, and is guaranteed to recover a solution of the actual Einstein 
equation R^i, = in coordinates defined by ^^ = 0. 

In the Riemannian elliptic context we are interested in the situation is 
a little more complicated. As we shall see later we must supply data on 
any boundaries, give an initial guess and essentially hope for the best. Our 
initial guess and our reference metric will typically be far from a Ricci flat 
solution. Hence we cannot consider the Harmonic Einstein equation only in 
the neighbourhood of a Ricci flat solution, as the entire problem is to flnd 
that solution. One might then imagine that the situation is hopeless in the 
elliptic context, that solving i?^ = is a quite unrelated problem to solving 
Rfiu = 0. This is not the case as we now discuss. 



2.2 Ricci flat solutions and Ricci solitons 

A solution to the Harmonic Einstein equation R^u = ^ laiv) with non- 
vanishing ^ is called a Ricci soliton. Obviously we are interested in Ricci 
flat solutions rather than solitons. As we proceed to discuss, fortunately the 

precisely the freedom to choose fixed charts on the manifold, and having made this choice 
there is no gauge freedom left. If the manifold has a boundary (or fictitious boundary as we 
discuss later) one may freely use a chart adapted to this boundary, so that the boundary is at 
some constant coordinate location, without being concerned that the boundary position in the 
harmonic coordinates might be something one must solve for as in some other gauge choices. 
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existence of solitons is rather constrained provided we choose our boundary 
conditions appropriately. 

Suppose we have boundaries or asymptotic regions in our problem and 
prescribe some data for the metric there, compatible with the ellipticity of 
the Harmonic Einstein equation. This data for the metric defines certain 
behaviour for the vector field ^. Consider as an example a manifold with 
a boundary. Taking local coordinates near the boundary we may write the 
metric as, 

ds'^ = a^dw'^ + -iij{dx' + /3'dw){dx^ + p^dw) (2.11) 

where the boundary is located at w = 0. Geometrically we might expect 
to try to fix some boundary condition involving the induced metric 7jj|^=o 
and the extrinsic curvature, Kij = ^ {d^Jij — 2V(j/3j)) |^=o (where the 
covariant derivatives and metric contractions are taken with respect to 7). 
We may regard 7ij as the geometric Dirichlet data, and Kij which involves 
the normal derivative of 7jj as the Neumann data. However 7jj and Kij 
only have D{D — l)/2 components, whereas the full metric g^u has D{D + 
l)/2. The Harmonic Einstein equation requires elliptic data for all these 
components of (7^1/. How should we fix the remaining D conditions? Since 
we are interested in Ricci fiat solutions of i?^ = where ^'^ = we had 
better enforce this. The additional conditions that ^^ = precisely give the 
extra D conditions for the metric componentsjj 

For the black holes discussed in this book we are not concerned with man- 
ifolds with boundary, but rather with asymptotic conditions. As we shall see 
explicitly later, we may impose an asymptotically flat or Kaluza-Klein con- 
dition also so as to ensure that ^ — )• at infinity. However this boundary 
example nicely illustrates that in order to give data for the Harmonic Ein- 
stein equation one should ensure that the data is consistent not only with 
imposing the geometric condition of interest, but also ensuring that ^ van- 
ishes. Obviously if one took boundary conditions that did not allow S^ to 
vanish then generally one might still be able to solve i?^ = but should 
not expect to find Ricci fiat solutions but only solitons. 

Consider a metric g and boundary/asymptotic conditions that we impose 
on it consistent with ellipticity. These boundary conditions impose certain 
behaviours on ^ - for example ^ = at a boundary as in the example above, 
or ^ — )• in an asymptotically fiat or Kaluza-Klein region. Then consider 

^ Very interestingly Anderson has shown that one cannot impose ^ = and require Dirichlet or 
Neumann boundary conditions (i.e. fixed induced metric or extrinsic curvature) as these are 
not well posed boundary conditions. Instead one may fix the conformal class of 7ij and the 

: 0.1231 
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the linear vector operator, 

v;^v^d; + R; (2.12) 



on this background so that we may write the Bianchi identity (2.10) as 
P • ^ = 0. Take a vector field x which has the same boundary or asymptotic 
behaviour as .^ - for example, x = ^ on the boundary described above, or 
X — )• in an asymptotically flat/Kaluza-Klein region. In order to find Ricci 
flat solutions to the Harmonic Einstein equation one should ensure that the 
metric boundary conditions lead to conditions on x such that the linear 
elliptic vector problem 2? • x = is well posed and admits the trivial solution 

x = o. 

Take the example of an asymptotically flat or Kaluza-Klein black hole. In 
the Riemannian signature with periodic time there are no boundaries and 
only the asymptotic region. The problem D-x = with the condition x ~^ 
asymptotically is indeed well posed and x = is obviously a solution. 

Note that we do not necessarily require the boundary condition that ^ ^- 
on all boundaries. For example, at a 'fictitious boundary' (see later section 



2.5) there are conditions as we shall see later where the tangential compo- 
nents of ^ have vanishing normal derivative (a Neumann condition), rather 
than being forced to vanish. This still gives a well posed problem P • x = 
consistent with a trivial solution. 

Now that we have boundary conditions consistent with obtaining a Ricci 
flat solution, we may use the Bianchi identity to control the existence of 
solitons. A necessary condition for a soliton to exist is that 2? • x = must 
admit a non trivial solution for x- Alternatively we may say that the vector 
operator D with the boundary conditions given by the behaviour of S, must 
have a non-trivial kernel. Now the kernel of such an operator should be finite 
dimensional, and may certainly be trivial. The necessary condition that D 
have a non-trivial kernel highly constrains the possibility of the existence 
of Ricci solitons. Indeed long ago Bourguignon [26] showed that there are 
no Ricci solitons on a compact manifold without boundary for any choice 
of vector ^. In certain cases such as for asymptotically flat or Kaluza-Klein 
metrics, as we shall discuss later we may prove the kernel is again trivial and 
solitons cannot exist for the DeTurck choice of ^ we are interested in. One 
arrives at the surprising conclusion that despite the fact that one is solving 
the Harmonic Einstein equation, which naively is quite different from the 
Einstein equation, in certain situations of interest the only solutions are in 
fact Ricci flat ones with the gauge condition ^ = imposed. 

Even if solitons do exist this is not a problem in principle. Since the Har- 
monic Einstein equation is elliptic, for well posed boundary conditions on the 
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metric we expect that solutions are locally unique. Hence a solution cannot 
be continuously deformed into another solution without suitably adjusting 
the boundary conditions. If there exists a Ricci flat solution, then any Ricci 
soliton solutions cannot be arbitrarily 'nearby' to it. Hence numerically it 
should always be possible to distinguish the solution of interest from the 
solitons. An obvious test is to compute the vector field (, and see if it is zero. 
In particular let us define the scalar = S,^£,n so that it gives the norm of 
^. For a Riemannian manifold vanishing ^ is necessary for the function (p to 
vanish. Hence we may check the magnitude of (j) computed for our solution, 
and if it is anywhere non-zero then the solution is a soliton and we should try 
to find another one. An important question in practice is how many soliton 
solutions there are. If one were trying to find a single Ricci flat solution in 
a vast forest of solitons then such an approach may be impractical. If it is 
a single Ricci flat solution in a small spinney it is practical. Obviously if no 
solitons exist then it is ideal. 



2.3 Asymptotically flat or Kaluza-Klein solutions 

Treating our static black hole problem as a Riemannian boundary value 
problem we must impose boundary conditions asymptotically. Our discus- 
sion now follows that developed in [19J with Figueras and Lucietti. We im- 
pose that the Riemannian geometry has an asymptotic region such that the 
metric approaches the direct product Sl x A4, where Sl is the Euclidean 
time circle with length f3 (corresponding to temperature T = 1//3) and Ai 
is a Ricci flat manifold, the obvious choices being Euclidean space M^"^ or 
•^D-2 ^ ^1^ ^]^g product of Euclidean space with a circle of length L. The 
former choice is appropriate for the Euclidean continuation of asymptoti- 
cally flat Lorentzian spacetime, and the latter would be for asymptotically 
Kaluza-Klein spacetime. 

Let us consider the case where we require the black hole to be asymptot- 
ically flat. We require the manifold to have an asymptotic region where for 
some large R the metric behaves as, 

ds'^ = g^ydx^dx" = dr^ + 6ijdx^dx^ + 0{r''P) 

digf^u = 0{r-P-') , didjg^, = 0{r-P-^) (2.13) 

for all r > i? and for some positive p. Here x* are the usual Euclidean 
coordinates with r = ^bijX'^xK For a Ricci fiat solution we expect p = D — 3. 



We require that the reference metric g which defines the vector ^ as in (2.9) 
is also asymptotically flat, so that in the same part of the manifold this 
reference metric also behaves as above. 



2 Static vacuum black holes 13 

In practice it is convenient to compactify the radial coordinate in this 
asymptotic region. Taking p = 1/r, we write the metric as, 



„2 1 

ds'^ = N^dr'^ + — dp2 + —habidO" + w''dp){de^ + w^dp) (2.14) 



" ^ iV 0,7 -|- 

p- ■ P' 



and require that, A^ = 1, a = 1, w^ = and hah = ^ab at p = 0, where 
^ab is the unit {D — 2)-sphere metric. These appear as Dirichlet boundary 
conditions at p = for the components of g, although we should recall 
that p = is really a regular singular point of the PDEs, since this is an 
asymptotic region rather than a boundary at finite distance. 

One may compute straightforwardly that taking an asymptotically flat 
reference metric, then, 

^- = 0{r-P-^) e = 0{r-P-^) . (2.15) 

We see that the norm (p = £,^^^ ~ 0{r^^^^^) so that the vector ^ goes to 
zero length asymptotically for any positive value of p. Thus these asymptotic 
boundary conditions are consistent with the linear elliptic problem D-x = 
discussed in the previous section being well posed and having trivial solution. 
The asymptotically Kaluza-Klein case proceeds in exactly the same manner, 
where one now requires that the metric asymptotes to, 

g = dT'^ + Sijdx'dx^ + dy"^ + Oir-'P) 
dig^u = 0{r~P-') , did.gf,, = 0{r-P-^) (2.16) 

for some positive p where y is the compact Kaluza-Klein circle, and has 
period L. One obtains the same behaviour, (j) ~ 0(r~^^~^), as in the asymp- 
totically flat case. 



2.4 A maximum principle and the non-existence of solitons 

We have discussed above that a static asymptotically flat or Kaluza-Klein 
black hole may be thought of as a smooth Riemannian manifold with no 
boundary except for the asymptotic region. We have shown that the vector 
^ — )• in this asymptotic region provided the reference metric shares the 
same asymptotics. We now turn to the important question raised earlier of 
whether we can control the existence of Ricci solitons with such boundary 
conditions. Following work with Figueras and Lucietti [19] we will now show 
that solitons cannot exist as we claimed earlier. This implies that solving the 
Harmonic Einstein equation is equivalent to solving the Einstein equation 
together with imposing the generalised harmonic gauge condition ^ = 0. 
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Contracting the Bianchi identity (2.10) with the vector and using the 
sohton equation yields, 

V'^cP + edf,<P = (V^e.)(V^r) > (2.17) 

where again (p = ^^S,f_i is the norm of the vector field and we have used 
the fact that for a Riemannian geometry the right-hand side is non-negative 
everywhere. Suppose now we have a solution to the Harmonic Einstein equa- 
tion. Consider a function / obeying the linear elliptic equation, 

VV + e''V>0- (2-18) 

Now this equation enjoys a maximum principle, which states that if / is 
non-constant, it must attain its maximum value on the boundary of the 
manifold. |j Furthermore, if a maximum exists at the boundary the outer 
normal derivative of / at this maximum is strictly positive. 

Since the function (f) > 0, and since a Ricci soliton must have (j) ^ 0, 
then a necessary condition for a solution to be a soliton is that (p is either 
constant and non-zero or it has a maximum somewhere. If (/> is constant 
and non-zero, then this implies V^^^/ = and so the vector ^ is covariantly 
constant, and in fact whilst the solution is a soliton it is also Ricci flat. Thus 
a challenge to Ricci flatness means that (p niust be non-constant and hence 



have a maximum somewhere. Since replacing / with (p in (2.18) we obtain 



the condition (2.17), and for (2.18) we have a maximum principle, we see 
that if (p is non-constant then its maximum must be at the boundary of the 
manifold, with positive normal outer gradient, or in an asymptotic region. 

In the case of no boundary, we recover the result that a solution must 
be Ricci flat. For black hole solutions which necessarily have some form of 
boundary, either at flnite proper distance from the horizon or asymptotically 
at infinite proper distance, we see that the existence of a soliton is intimately 
tied to the precise form of the boundary conditions imposed there. 

Let us consider the asymptotically flat or Kaluza-Klein boundary condi- 
tions discussed in the previous section. Assuming we flnd a solution to the 
Harmonic Einstein equation compatible with those asymptotics we calcu- 
lated that (p = l^p = 0(r~^P~^) and hence for the very weak requirement 
that p > 0, this implies (/>—)• asymptotically. Provided there are no other 
boundaries, a simple application of the maximum principle rules out the 
existence of Ricci solitons as we claimed earlier£J 

^ See for example 27: , or for a statement of results specifically on Riemannian manifolds see 

m- 

"^ Assume for contradiction that there is a Ricci soliton so that 4> j^ 0. Consider some R' > R, 
so that on the closed surface r = R' in the asymptotic region we have (j) < C jR!"^^^^ for some 
constant C. Then everywhere in the interior of this surface the maximum principle states that 
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2.5 Isometries and 'fictitious' boundaries 

Let us now turn to a more practical issue, namely using isometries to reduce 
the effective dimension of the PDEs. A PDE problem in high dimensions 
rapidly becomes intractable. Thus if we have isometries in our geometry 
it is important to manifest them explicitly using adapted coordinates and 
hence benefit from the storage savings of not having represented the met- 
ric in the isometry directions. The example of the 5D Kaluza-Klein black 
holes discussed in Chapter 4 is a case in point. This 5D problem is reduced 
to 2 effective dimensions when one takes into account the static ^(1) and 
rotational SO{3) isometries. 

Suppose we wish to numerically represent a function that has spherical 
symmetry in Euclidean space M™. Using Cartesian coordinates, so ds^ = 
Sijdx^dx^ , the origin of the symmetry x* = is not a special point and 
the function is simply smooth there. If we use coordinates adapted to the 
symmetry, spherical polar coordinates ds'^ = dr'^ +r'^dQ,'^, the function only 
depends on the radial coordinate r giving the saving in storage. However, 
we find that the polar coordinate chart breaks down at the origin, and 
we must treat the point at the origin effectively as a boundary point - we 
term this a 'fictitious boundary'. Of course we may deduce the required 
boundary conditions in the polar chart simply from the requirement that 
in the Cartesian coordinates the function is smooth, so that /(x) is a C°° 
function. Then using the fact that r^ = dtjx'^x^ this implies a function that 
only has radial dependence must be a smooth function in r^, i.e. / = /(r^) 
isC°^. 

There are two important fictitious boundaries that typically arise when 
considering static black holes. Firstly associated to the Euclidean static ^(1) 
isometry we have that the horizon is precisely the place where the time circle 
vanishes, and hence the isometry has a fixed action. In this case we may write 
the metric in polar coordinates as, 

ds^ = Adr^ + r^Bdr^ + rCadrdx" + hahdx^dx^ (2.19) 

where the component functions only depend on the radial and x coordinates, 
and are independent of the time direction r which we recall has period 
T ^ T + Itx I K. Using the Cartesian coordinates X and Y mentioned above, 
so X = r cos KT and Y = r sin kt, then one can show the metric components 

4> < C / R''^^'^^ . Under the assumption that <f> j^ 0, we must be able to find a point p where <j> 
is non-zero, and denote its value <po, so </>o = ipip)- Choosing R' large enough, we may always 
ensure that the point p is in the interior of the portion of the manifold bounded by the 
surface r = R' and also that <j>o > C/R''^^^^ . Since we know in this interior 4> < C/R'"^^^^ we 
are lead to a contradiction and must conclude that </> = everywhere as a consequence of the 
maximum principle. 



16 Numerical construction of static and stationary black holes 

in these coordinates are smooth functions of X and Y at the fixed point 
r = provided that A, B,Ca,hab are smooth functions of r^ = X"^ + Y'^ 
and the coordinates x, and furthermore that k?A = B at r = for some 
constant k > 0. Here the constant k gives the surface gravity of the horizon 
with respect to d/dr. 

Since there is not a real boundary at the horizon our maximum principle 
must ensure no maximum of (p can occur there. Treating the horizon in 
static adapted coordinates we see the horizon as a fictitious boundary. We 
may simply compute that ^'' = dr^"" = at r = 0. Now ^"^ = everywhere 
and hence dr(j)\r=o = 0. Recall the maximum principle implies a maximum 
at a boundary requires (p has positive outer normal gradient. We see since 
dr(f) = at r = 0, the maximum principle does indeed rule out a maximum 
at this fictitious boundary. 

Secondly one typically has rotational axisymmetry in the problem and an 
associated SO{n) isometry group. The axis of this symmetry is then the set 
of points fixed under the isometry. In this case we write the polar metric as, 

ds'^ = Adr'^ + r'^Bdn'^ + rCadrdx" + habdx^dx^ (2.20) 

where dVt^ is the line element on a unit (n — l)-sphere. Transforming to 
appropriate Cartesian coordinates one finds a smooth metric again provided 
A^ B, Ca, hab are smooth functions of r^ and A = B at r = 0. Note that 
there is no free constant as in the case of the vanishing U{1) as we have 
already chosen the sphere dQ'^ to have unit radius. As above, no maximum 
can reside here as dr(f> = at r = 0. 



2. 6 Solving the Harmonic Einstein equation I: Ricci flow as 

local relaxation 

We shall now consider the canonical method to solve an elliptic system, and 
find that in fact for the Harmonic Einstein equation this can be thought of 
in the continuum as the famous Ricci flow. 

Suppose we consider the Laplace equation, V'^ip = 0, in some flnite region 
U of D-dimensional Euclidean space, and wish to solve this boundary value 
problem with some appropriate boundary conditions (e.g. Dirichlet). The 
simplest numerical approach is to represent the function ip using real space 
flnite difference and use local relaxation to flnd a solution. Let us take the 
canonical coordinates (xi, X2, . . . , xd) on R-^, and consider the rectangular 
lattice of points with lattice spacing A at positions (mi A, m2A, . . . , m/jA) 
for integers mi,m2 ■ ■ ■ , mjj. Denote the set of points {pi} in this lattice that 
lie in the interior of the domain U as L, where i = 1, . . . ,N labels these 
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points. Real space finite difference represents tlie function ip by storing tlie 
set of values {ipi} where ipi = '4){p,i). 

The simplest way to represent the Laplace equation is by second order 
finite difference. Consider a point pi £ L. We denote the 2D nearest neigh- 
bour points in this rectangular lattice as Pi->j, where j = 1,2,..., 2D. Then 
we may approximate, 

V'V'U ^ ^ I -(2I))^, + ^V.->, j • (2.21) 

By making A smaller, we approximate this Laplacian increasingly well, and 
we now have a finite representation of the continuum Laplace equation. Note 
that for points in L that neighbour points in the boundary or exterior to U, 
then we consider the values of those neighbours to be fixed by the boundary 
conditions. 

We then proceed by solving this finite problem. The classic method is 
that due to Jacobi, a local iterative procedure also known as relaxation. 
Let us imagine a sequence of guesses to the solution, {il^^ } for integer 
^ = 0, 1, 2, . . .. The Jacobi method states that given a guess {ipl } we can 
improve it by computing a new guess {t/j- } as 



^ (I*"') 



-/■r" = ^ I > >;-ii I (2.22) 



The idea is that one takes an initial guess {ipl } satisfying the boundary 
conditions, and then iterates Jacobi's improvement method. For the simple 
Laplace equation one will reach a fixed point of this iteration which is a 
solution of the finite Laplace equation above. 
We may rearrange the equation above as, 

2D 

Y.4%-imi^r'\ (2.23) 

but now this may be viewed as the finite differencing of the diffusion equa- 
tion, where now ■i/' is a function of x and a fiow time A, and, 

^^ = VV(A,.) (2.24) 

where we discretize time similarly to space, so that tp^ = ^(A(5)|p- where 
(5 = ^ gives the discretization spacing in flow time. The right-hand side is 




^(#"^^-#0- (^-2^) 
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approximated as above, and the left-hand side is approximated using forward 
Euler differencing, 

9^ X=A5 

Thus we see that the simplest local relaxation method used to solve the 
Laplace equation is in fact simply diffusion on spatial and temporal scales 
must larger than the lattice scale A. Instead of thinking about Jacobi local 
relaxation and finite difference of the elliptic problem, we can describe this 
process more formally as the parabolic continuum diffusion problem. Then 
we take an initial guess, and act on it with diffusion until we reach a fixed 
point of the diffusion flow, and this will solve the original elliptic problem. 
We may use the Jacobi method for the more complicated elliptic Harmonic 
Einstein equation. Instead of a single function ^, in a chart there are all the 
components of the metric tensor g^^, to solve for. We may represent the 
metric in a chart by again taking the same lattice L in M , and we store the 



set of values {gp,u)i = 9fj.u{Pi)- Recall that from (2.8) the Harmonic Einstein 
equation has second derivative structure, 

<. = -Ig'^'^dad^gf^u + T^u (2.26) 

where T^i, represents the terms with lower numbers of derivatives. Jacobi's 
method is then the discretisation of the continuum equation, 

^^ = -2R^. (2.27) 

where the right-hand side is approximated using second order finite differ- 
ence, and the left-hand side is differenced using the forward Euler method. 
We may think of Jacobi as approximating the Harmonic Einstein equation 
as a set of Poisson equations, where the sources in the Poisson equations are 
the lower derivative terms T^u- It attempts to solve the Poisson equations at 
each point in isolation by considering the neighbours and sources as fixed. 

In fact there are many variants of the classical Jacobi method of relax- 
ation for elliptic systems. For example Gauss-Seidel speeds up Jacobi by 
updating the values of the points in situ, rather than storing a whole new 
approximation ip('^~^^) given a previous one Tp^ '. This is again another dis- 
cretisation of the continuum diffusion problem although with an exotic time 
derivative differencing. Over and under relaxation schemes simply adjust 
the diffusion constant. Multi-grid methods work by concurrently solving the 
equation on multiple lattices, cascading information from low to high resolu- 
tion and back, in order to allow information to propagate quickly over long 
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distances. These methods give great enhancement of speed, but on scales 
larger than the lowest resolution lattice, again act as continuum diffusion. 

Thus whichever local relaxation method one chooses, one may consider 
the problem to be continuum diffusion on large enough scales, and for the 
Harmonic Einstein equation, the canonical diffusion is the flow in time A 



given in equation (2.27), which explicitly is 



dX 



2R^, + 2V(^e.) ■ (2.28) 



In fact this is precisely the Ricci-DeTurck flow, where the second term is 
an infinitesimal diffeomorphism generated by the vector field ^, so that the 
flow is diffeomorphic |j to the Ricci flow 

^^ = -2R,u (2.29) 

introduced by Hamilton as a tool in geometric analysis, and which has gained 
fame for its role in proving the Poincare conjecture. For an introduction to 
Ricci flow see for example ^29j. DeTurck proved that Ricci flow was a well 
posed parabolic flow by realising that one can add the diffeomorphism term 



and with the choice in (2.7) explicitly render the flow equation parabolic. 

Again on scales larger than the lattice resolution used in the local re- 
laxation, we may view the methods discussed above formally as the Ricci- 
DeTurck flow. We may regard this parabolic flow as a continuum algorithm 
to solve the Harmonic Einstein equation. We give some initial guess for 
the parabolic flow. A fixed point of the fiow is a solution of the Harmonic 
Einstein problem. Hence, by simulating the flow for sufficient flow time, we 
might hope to approach a fixed point as closely as we require. 

One beautiful consequence of this is that whilst we required some choice of 
reference metric to define the vector field ^ and render the Harmonic Einstein 
equation elliptic, and hence well posed as a boundary value problem, in 
fact the Ricci-DeTurck flow is diffeomorphic to Ricci flow which makes no 
reference to ^. Thus given some initial guess metric, whilst different choices 
of reference metric will change the path taken in the space of metrics by the 
Ricci-DeTurck flow, the path taken in the space of geometries (i.e. metrics 
modulo diffeomorphisms) is always the same. 

Provided one chooses the reference metric to share the same isometries 
as the metric, for example the static symmetry, then the Harmonic Einstein 

* In the presence of boundaries we should ensure that the normal component of ^ at a 
boundary vanishes in order that these two flows are diffeomorphic. This ensures that the 
diffeomorphisms generated by ^ along the flow act to preserve the boundary points. 
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tensor is also symmetric under these isometries. This imphes that the Ricci- 
DeTurck flow preserves the isometries of the metric. 

Ricci flow has very nice properties geometricahy. EssentiaUy it is diffusion 
for geometry, and locahy tries to smooth out curvature. Suppose we have a 
Ricci flat solution, ^r^j^, and wish to consider the Ricci flow of a perturbation 
to this, h^i^. Then from (2.4) we see, 



dh 



flU 



dX 



-2Al V - 2V(^i;^) , (2.30) 



where Al is the Lichnerowicz operator on the background g, and the last 
term is simply an inflnitesimal diffeomorphism generated by v. For Euclidean 
space, this flow is then diffeomorphic to a flow where each component of the 
metric simply diffuses as, -qj^ = ^°^^dadph^v. Hence we see that Euclidean 
space is stable to linear perturbations. However beyond this diffusive be- 
haviour of linear perturbations it has very interesting non-linear properties. 
For example, it wishes to collapse regions of positive curvature as can be 
seen by the Ricci flow of a round sphere, whose radius shrinks linearly in 
flow time, reaching a zero size in flnite time. Of relevance for us, it has been 
rigourously proven that Ricci flow exists and preserves asymptotic flatness 
for short times [30] . 

An important property of the Ricci flow is that we see from above a 
Ricci flat solution is a stable flxed point only if the operator A^ is positive. 
We may assume there are no zero modes, as these should be flxed appropri- 
ately by boundary conditions, which should ensure a locally unique solution. 
However A^ may not be positive in general, so that there exist one or more 
eigenfunctions with negative eigenvalue. At late times such eigenmodes ex- 
ponentially grow in time so that the perturbation will flow one away from 
the flxed point in these directions. 

A fundamental property of static vacuum black holes is that the positivity 
of Ax, for their Euclidean continuation is related to their thermodynamic 
behaviour, since the Euclidean action is simply related to their free energy. In 
the asymptotically flat case this has been made precise in [3l] . In many cases 
static black holes of interest (for example, all those in Kaluza-Klein theory 
discussed in Chapter 4) possess negative modes of A^, the canonical example 
being that of the asymptotically flat Euclidean Schwarzschild solution which 
has the single negative mode discovered by Gross, Perry and Yaffe [32j. 

We have seen that standard local relaxation methods can be viewed as 
Ricci-DeTurck flow on large scales, which is diffeomorphic to Ricci flow. 
Thus for all these methods, one cannot simply start with an initial guess 
and flow to the black hole flxed point solution if it possesses negative modes. 
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Starting with an initial guess close to the solution one will flow towards the 
fixed point in nearly all directions in the space of perturbations of the fixed 
point, but will generically then veer off along the direction(s) tangent to the 
negative mode(s). 

Since many solutions of interest have negative modes one might imagine 
that local relaxation hopelessly fails to provide an algorithm to find solutions 
of the elliptic Harmonic Einstein problem. This conclusion however is too 
quick, and in principle Ricci flow and hence local relaxation may still be used 
but the method must be modified slightly. Suppose the fixed point we are 
interested in has a single negative mode of Al and up to diffeomorphisms 
has no zero modes so that the fixed point is locally unique. Examples of this 
include localised black holes in Kaluza-Klein theory discussed in Chapter 4. 
Let us call the fixed point go. Then locally about go the space of geometries 
(meaning metrics up to diffeomorphisms) is infinite dimensional. There are 
two special fiows that emanate from go along the negative mode direction. 
Let us denote the negative mode perturbation h at go. Then these two flows 
are g±{X) such that g± c^ go i: e^"^^ '^/i as A — )■ — oo, where A^^/i = —v^h so 
that v^ is the magnitude of the negative eigenvalue. The perturbation /i is a 
tangent vector to the space of geometries at ^o- A basis for the tangent space 
at go is given by h together with the positive eigenmodes of A^. The positive 
eigenmodes are tangent to a codimension one surface, E, which contains ^fo 
and is closed under Ricci flow. Starting from any point in S near to go one 
remains within this surface under the action of the flow, and will flow to the 
flxed point reaching it asymptotically. 

The problem is then to generate an initial guess contained in E since this 
will flow to the fixed point. Consider a one parameter family of geometries 
gia), where a is the parameter. It is generic that the curve gia) will intersect 
E. Suppose this occurs at a = a^,. Then for a> a^, but close to a^ one will 
initially flow towards ^o and then be carried away in the direction of the 
negative mode. In the limit that a — )• a*, the flow will at late times follow 
that of either g^ or g^. Let us assume it follows 5+ for a > a-),. Then 
conversely for a < a* one will approach go and then deviate away in the 
opposite sense, flowing away close to the flow g-. 

Thus we see that for a close to a* there is a critical behaviour associated 
to the unstable fixed point go. By scanning the values of a one can hope to 
see this critical behaviour, and if one can identify whether one has fiowed 
in the (7+ or </_ direction, one can simply automate a tuning of a to get as 
close to a^ as required. Then one has flows that get very close to go for a 
long period of flow time before flnally succumbing to the negative mode and 
flowing away. In principle one can get as close to the flxed point as desired. 
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For a Schwarzschild black hole the flows g± generated by the negative 
mode either expand (say g^) or contract (g_) the horizon. The flow g- has 
been shown to continue to shrink the horizon to a finite time singularity, 
whilst under the flow g^ the horizon grows without stopping [33j. Thus 
given a flow it is very simple to see which side of S that flow is on, and 
hence tune the parameter a to reach a>,. Similar behaviour is seen for the 
Kaluza-Klein localised black holes, where one flow pinches the horizon to 
zero size and a flnite time singularity, and the other expands it until its 
poles touch and again a singularity is reached jH| . 

In principle this method may be extended to a case with A^ negative 
modes. Then an N parameter family of initial data must be tuned in order 
to reach the flxed point. In the codimension one case, so that S partitions 
the space of geometries locally about go, it is easy to see which 'side' of S 
one starts on. This is not true in the higher codimension case, and one must 
simply search the space of parameters until one locates the critical point a*. 

Let us summarise this discussion. Local relaxation is the simplest method 
to solve elliptic PDEs, and can be applied to the Harmonic Einstein equation. 
On large scales we may think of relaxation from a continuum perspective 
as Ricci-DeTurck flow, which is diffeomorphic to Ricci flow. Hence this ap- 
proach has the beautiful geometric property that the trajectory taken by 
the flow is independent of the reference metric and hence the gauge flxing. 
However many black holes are unstable flxed points Ricci flow, and for these 
relaxation or Ricci flow may still be used to flnd these solutions, but for a 
solution with N negative modes of its Lichnerowicz operator, one must find 
a suitable A'^ parameter set of initial data, and tune these A^ parameters 
in order to fiow or relax to a solution. We emphasize that there do exist 
interesting solutions with Killing horizons which are stable under Ricci fiow, 
such as those in AdS/CFT where the boundary metric is a black hole |19] . 



2.7 Solving the Harmonic Einstein equation II: Newton's method 

We have seen that whilst the simplest relaxation methods to solve elliptic 
systems have an elegant geometric behaviour on large scales, it is difficult 
to find many black holes of interest which are unstable fixed points of these 
methods due to having Euclidean negative modes. With one such negative 
mode these methods are still practical. For more they become increasingly 
hard to use. 

Fortunately there is a second standard technique to solve these elliptic 
systems, namely Newton's method (also known as the Newton-Raphson 
method). As we shall see this approach is considerably more complicated 
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to implement, and lacks the geometric elegance of relaxation, so that the 
behaviour of Newton's method will explicitly depend on the choice of ref- 
erence metric. However, the advantage of Newton's method is that it is 
insensitive to the stability of the fixed point. In fact the basin of attraction 
of Newton's method can be rather small in practice, and thus a combina- 
tion of Ricci flow or relaxation to get close to the fixed point, followed by 
Newton's method to hone in on it can be the best strategy. 

Unlike relaxation, Newton's method is inherently non-local. Let us again 
imagine discretising our system using finite difference as above. At each 
lattice point in a chart we will have the various components of the metric. 
Globally there will be a finite set of numbers {gAi} that will give the finite 
difference approximation to the metric g^u{x), where the index M includes 
both the lattice point in a given chart and the component of the tensor. 
Likewise we may represent the Harmonic Einstein tensor in the interior of 
the manifold with the same index structure. Then the Harmonic Einstein 
equation is given as R^i{g) = 0, which can be thought of as a finite set 
of coupled non- linear equations in the variables gM- The canonical way to 
solve such a system is by the multidimensional generalisation of Newton's 
method. 

If we perturb the metric g as g + eSg, the Harmonic Einstein tensor goes 
as, 

R^ig + e5g) = Rf,{g) + e 0{g)^^6gN + 0{e^) , (2.31) 

where the matrix 0{g)jff is the linearisation of Rfj- Begin with an initial 
guess gj^ . Then Newton's method iteratively improves a trial metric g)^' 
as,E] 

gir'^=g^j^^-m9^^YX^RM^^). (2.32) 

As with the one dimensional Newton method this moves along the tangent 
of the equations to find a solution. Near a solution it will very quickly con- 
verge to that solution. However, the basin of attraction may be rather small 
in practice, and outside of this iterations of Newton's method will usually 
diverge and give singularities. 

As for the Ricci flow method, provided the reference metric is chosen to 
have the same isometrics as the metric, then the Harmonic Einstein tensor 
will be symmetric under these and the Newton method will act to preserve 
these isometrics. 

^ It is sometimes useful to take 'smaller' steps, with gy^+^) = g(^> — eO{g^ ')~^ ■ R^ {g^ ') 
for some e with < e < 1, particularly in the first iterations if the initial guess is not very 
close to the solution. 
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This method has the important advantage over the Ricci-DeTurck flow 
method that it is not sensitive to negative modes of the Lichnerowicz op- 
erator. However it does assume that the hnear problem O ■ V = R^ can 
be solved for the vector V. In practice robust methods exist to solve such 
(finite dimensional) linear systems, such as biconjugate gradient, which are 
insensitive to the spectrum of O, provided there are no zero modes which we 
assume for well posed boundary data. Thus a single initial guess will suffice, 
rather than having to tune a family of initial guesses. 

We see that the implementation of Newton's method is considerably more 
complicated than that of the relaxation/Ricci flow methods. Another impor- 
tant disadvantage of the Newton method over the Ricci-DeTurck flow is that 
it is not geometric in the sense that the path taken by the algorithm in the 
space of geometries will depend explicitly on the choice of reference metric. 
This implies that the basin of attraction of a solution, which in practice 
may be rather small, will also depend on this choice of reference metric. 
Sometimes it is actually convenient to use a combination of the Ricci flow 
method together with the Newton method. The Ricci flow method is rather 
robust and can quickly get one reasonably close to a fixed point. It is tun- 
ing the flows very close to the fixed point that becomes difficult and time 
consuming. However, once reasonably close, one can simply use the Newton 
method to quickly find the precise fixed point. 



2.8 An illustrative example 

In order to illustrate many of the points discussed above we will now give a 
very simple example, where we already know the answer. We will consider 
finding the 4D Schwarzschild solution using the various techniques above. 
We will assume spherical symmetry, and so the problem really is one in- 
volving only ODEs but we shall treat it in an identical manner to the much 
more complicated PDE problems we are really interested in. This example is 
simple enough that we can be very explicit about the implementation which 
we detail below. 



10 



We note that our earlier maximum principle argument states that in this 
asymptotically fiat case no soliton solutions should exist and hence any 
solution to the Harmonic Einstein equation must be Ricci fiat. We will cover 
the manifold with one chart, and use the fact that the solution is static and 
spherically symmetric, adapting coordinates to these symmetries. We will 

^^ We will make available a very simple Mathematica notebook that implements the relaxation 
/ Ricci flow and the Newton algorithms in this toy example. Hopefully this provides an entry 
point for those interested in thinking about the more complicated problems of interest. This 
will be found at; http://www3.imperial.ac.Uk/people/t .wiseman 
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choose a radial coordinate r, and choose the horizon to be located at r = 0, 
and infinity to be at r = 1. We continue to Euclidean time, and write the 
smooth Riemannian metric as, 

ds'^ = r^AdT^ + AfBdr^ + fCdn'^ , / = ^ K^^ (2.33) 

(1 — r^)^ 

where d^^ = dO"^ + sin 9d(fP'. Then A,B,C are functions of r, and the 
metric describes the general static spherically symmetric metric. Our choice 
of factors above means that for the topology of manifold we wish to describe, 
^, -B, C > on the domain r € [0, 1]. At r = we have a fictitious boundary 
as we have adapted coordinates to the static symmetry. From our previous 
discussion smoothness of the full Riemannian manifold implies that at r = 
then K^ = A/AB for surface gravity k, and A, B, C are smooth functions of 
r^. At infinity, r = 1, we fix A = S = C = 1. 

The Harmonic Einstein tensor is determined by the components i?^, i?^ 
and Rqq. We may discretize this system using finite difference by choosing 
(1+A^) lattice of points at locations rj = iA for i = 0, 1, . . . A^ with A = 1/A^, 
so that Ai = A{ri), and likewise for B and C. Our boundary conditions imply 
that An = Bjy = Cn = 1 . At the horizon we require smoothness in r^ which 
implies that we may deduce the boundary values Aq^Bq^ Cq in terms of the 
interior points. For small r the behaviour goes as constant plus quadratic in 
r, which gives ^o = (4Ai — A2) /3, and likewise for Bq and Cq. We choose 
to also impose the regularity condition A = B directly by taking Aq = Bq. 
This in turn determines Bi from the smoothness. In full we have, 

^0 = ^ (4^1 -A2), Bo = \ (4^1 - A2) , Co = ^ (4Ci - C2) , 

51 = ^1 + ^(^2-^2), An = Bn = Cn = 1 (2.34) 

and the vector qm = {Ai . . . ^Af-i, B2 ■ ■ ■ Bn-i-, Ci . . . Cjv_i} describes the 
metric subject to the above conditions. We finite difference derivative terms 
using simple second order differencing so that, 

drXi = —— {Xi^i — Xi^i) 

dlX, = ^ (Xi+i + X,_i - 2X{) (2.35) 

and then we may evaluate R^^. = Rjj^{r,i). 

Consider gjj to be an initial guess, and gj^^ to be subsequent iterations 
of improvement for A = 1,2, . . .. Then the Jacobi method, or equivalently 
the Ricci-DeTurck flow discretised in time using forward Euler differencing. 
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gives, 

rj^^'^ = rf^ -26R?A9^^\ , ^ = l,...,N-l (2.36) 

and similarly for B and C except that for B we have i = 2, . . . ,N — 1. The 
continuum Ricci flow time A for g^ > is then given as X = A6 and 6 = A^/2. 
For the Newton method one creates the (3A^ — 4) vector of equations, 

r>H — ( nH r,H nH nH nH nH \ (r, o7\ 

^M ~ l^TTl ■ ■ ■ ^TT AT-D ^rr2 ■ ■ ■ ^rrN-1^ ^901 ■ ■ ■ ^eON-lJ {■^■Ol ) 

which is a function of the (3A^ — 4) component vector qm- Then the lineari- 
sation Oj^ = dR^^/dgN is a square matrix, which can be inverted to solve 
the linear system required for the Newton method. 

Without loss of generality we choose k, = 1/2 by global scaling. Now we 
note at this point that A{r) = B{r) = C{r) = 1 is the Schwarzschild solution 
for n = 1/2. Whilst this is a toy example, it almost seems too trivial to find 
the solution in these coordinates, and hence to challenge ourselves we will 
choose the background metric such that for the Schwarzschild solution the 
metric functions are not simply constant! 

We know that Schwarzschild is unstable to Ricci flow with one negative 
mode. Thus let us choose a one parameter family of initial metrics, 

A = l-a{l-r^f , B = l-a{l- r'^f , C7 = ^ (l + r^) (2.38) 

parameterised by the constant a that satisfy our boundary conditions for 
K = 1/2. We choose the metric at zero flow time to be given by A, B and 
C above, and also take the fixed reference metric to be the same. We note 
that since C above is not constant then the actual Schwarzschild solution 
has non constant metric functions in the generalized harmonic coordinates 
the reference metric imposes. 

Under relaxation/Ricci flow we might hope there is a critical value of 
a, say a^ where we may approach the unstable Schwarzschild fixed point. 
Indeed one finds there is, with the critical value a^ ~ 0.72. In figure [L] we 
plot the size of the sphere at the horizon R = \/C\r=o against fiow time A for 
a number of flows with various a approaching a^ from above and below. We 
see that two very distinct behaviours are found for a > a* or a < a*, that 
are straightforward to identify, and hence tune to the fixed point. Needless 
to say, the fixed point itself is indeed Schwarzschild, and we see that for 
a ~ a^ the horizon does indeed tend to unit radius at late times as it should 
for Schwarzschild with k = 1/2. 

The Newton method efficiently hones in on the Schwarzschild solution 
provided one is in the basin of attraction of the fixed point. In fact for the 
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Figure 1. Figure depicting use of relaxation or the Ricci flow method 
to find the Schwarzschild solution. The plot shows the evolution of the 
radius of the horizon, R, as a function of flow time A for flows with a 
variety of values of the parameter a. The values of a are labelled next to 
the corresponding curve. After initial transient behaviour we see that for 
a < a^, ~ 0.72 the horizon shrinks (in fact to a singularity) at finite fiow 
time. Conversely for a > a^, the horizon grows indefinitely. Tuning close to 
a* one may approach the Schwarzschild solution as accurately as desired. 
We see the quantity plotted here does indeed tend to one (the Schwarzschild 
value for the chosen surface gravity) for a finely tuned flow. (The data 
presented was computed using a very modest A^ = 100; we note the value 
of a required to flne tune the data will depend on this discretization, giving 
a* only in the continuum limit). 

choice of initial metric and reference metric given above, a rather wide range 
of values of a all lie in the basin of attraction. For example, simply taking 
the initial guess with a = will quickly find the Schwarzschild solution after 
a handful of Newton iterations. 



3 Stationary vacuum solutions 

Generally we wish to be able to tackle stationary solutions and in this section 
we address how to extend the static methods discussed above to this case. 
We note that the classic 4D uniqueness theorems relied on formulating the 
stationary axisymmetric problem as an elliptic system [2j. Our task here is 
to formulate the general stationary vacuum problem as an elliptic system 
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using the covariant Harmonic Einstein equation approach, and then ensure 
that the Ricci flow and Newton method algorithms may still be applied. This 
will require us to tackle the problem from a manifestly Lorentzian point of 
view. This section is based on the recent work [34] together with Adam and 
Kitchen. 



3.1 Static solutions from a Lorentzian perspective 

Instead of immediately considering stationary spacetimes, it is instructive to 
first consider static spacetimes from a Lorentzian perspective. The Harmonic 
Einstein equation is not elliptic for a general Lorentzian manifold but rather 
it is hyperbolic, and without ellipticity one would not expect to be able 
to impose the various boundary conditions that physically we require in a 
well posed manner. However, consider a chart away from any horizon which 
manifests the static symmetry, 

ds"^ = g^^dxi'dx" = -N{xfdt^ + hij{x)dx'dx^ (3.1) 

so that A^^ > 0. With the choice that our reference metric is also static with 
respect to d/dt, so that, 

ds"^ = g^,udx''dx'' = -N{xfdt^ + hij{x)dx'dx^ (3.2) 

again with iV^ > and hij a smooth Euclidean metric, then i?^ shares 
the static symmetry. Then due to this static symmetry the Harmonic Ein- 
stein equation i?^ = thought of as PDEs for the metric components of g 
is invariant under an analytic continuation t ^ t = it. Hence we immedi- 
ately see that the Harmonic Einstein equation restricted to Lorentzian static 
metrics and reference metrics is elliptic. The behaviour of local relaxation 
or Ricci flow, and the Newton method will be precisely the same in either 
signature. 

Under a Euclidean continuation, with Euclidean time taken to be peri- 
odic, we could remove the boundary associated to a horizon. However, we 
noted that in practice one should take advantage of the static isometry and 
adapt coordinates to it, but then the horizon manifests itself as a flctitious 
boundary of such coordinates, analogous to the origin of polar coordinates. 
Boundary conditions at the horizon in these adapted coordinates are just 
derived from transforming to regular coordinates which do not manifest 
the isometry but do manifest smoothness. In Lorentzian signature we have 
no option but to think of the horizon as a boundary. However, since the 
Harmonic Einstein equations are independent of signature, the boundary 
conditions for a regular Lorentzian horizon are precisely the same as those 
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in Euclidean signature. Let us take coordinates in the base adapted to the 
horizon such that x* = (r, x") where r = is the horizon. Then we write, 

ds'^ = -r'^Vdt^ + Udr"^ + r Uadrdx"- + habdx'^dx^ (3.3) 

where the metric functions are functions of r and x". Changing to coordi- 
nates, 

a = rcoshKt, 6 = rsinhKt (3.4) 

provides a good chart covering the static Kilhng horizon, such that the met- 
ric components are smooth functions, provided that V, U, Ua, hab are smooth 
(C°°) functions of r^ and x'^, and, 

V = K^U (3.5) 

at the horizon r = 0, and k is again the surface gravity. The same conditions 
wiU apply to the reference metric which is also required to be smooth at r = 
0. Of course we just have exactly the same situation as in the Euclidean case 



in equation (2.19) in <:2.5 with time continued back to Lorentzian signature. 
In the Euclidean context it is clear that since the metric is smooth and 
without boundary at the horizon, i?^ must also be smooth there. The same 
is true in the Lorentzian case where the tensor i?^ shares the same regularity 
properties as the metric. This is simply seen by noting that in the coordinates 
(a, 6, x") the metric and reference metric components are smooth functions, 
and hence so are those of Rjj^- Transforming back to the static adapted 
coordinates (t, r, x") then gives, 

R" = -r^fdt'^ + gdr"^ + r gadrdx" + rabdx^dx^ (3.6) 

where f,g,ga and Tab are smooth in r'^,x", and in addition / = K^g. Thus 
in the Lorentzian picture we have the nice property that Ricci flow and 
the Newton method will preserve the regularity of the horizon boundary. 
Furthermore it will preserve the surface gravity of the horizon. Thus we must 
now view the horizon as a boundary and we are naturally lead to impose 
physical data there, namely the surface gravity with respect to d/dt. 



3.2 Stationary spacetimes with globally timelike Killing vector 

We begin our discussion of stationary spacetimes by considering the case 
of spacetimes with globally timelike Killing vector, and will argue that the 
Harmonic Einstein equation is elliptic. Of course we are ultimately interested 
in black hole spacetimes which violate such a condition, with the stationary 
Killing vector becoming either null on the horizon, or outside the horizon 
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at the boundary of an ergoregion. In the following section ^3.3 we consider 
more general stationary spacetimes which allow horizons and ergoregions. 

Consider the most general stationary metric with Killing vector T = d/dt, 
which we may write using coordinates adapted to the stationary isometry 
as, 

ds^ = -N{x) {dt + Ai{x)dx'f + hij{x)dx'dx^ . (3.7) 

Now under our assumption that T is globally timelike we have A^ > and we 
further assume that the function N is bounded. Physically this implies our 
spacetime has no Killing horizons, and also no ergoregions. Since del g^i, = 
—N det hij we see that provided the metric g is Lorentzian and smooth, so 
that detg^u < and bounded, this implies that det hij > 0. We may then 
regard this metric as a smooth fibration of time over a base manifold Ai 
so that {Ai, h) is a smooth Riemannian manifold with Euclidean signature 
metric hij. It is worth noting that this metric is not of the ADM form, 
but rather takes the form of a Kaluza-Klein reduction ansatz with respect 
to time. Thus hij does not give the metric of a constant time slice of the 
Lorentzian geometry. 

The second order derivative terms acting on the metric components g'^^/ 
in the stationary Harmonic Einstein equation go as, 

R% = -\9''^dadpg^, + ... = -hi'^didjg^, + ... (3.8) 

where . . . are lower derivative terms. We see whilst the metric g^u is indeed 
Lorentzian, since there is no dependence on the coordinate t, it is actually the 
metric hij that controls the character. We note the Kaluza-Klein form above 
ensures that the inverse metric g^^ in the base directions is simply given in 
terms of the inverse of hij. This immediately implies that the Harmonic 
Einstein equation i?^ = is elliptic since h is smooth and of Euclidean 
signature. 

We also require that i?^ is a tensor that is symmetric with respect to 
the stationary isometry T. Without this, Ricci-DeTurck flow and Newton's 



method will not consistently truncate to the class of stationary metrics (3.7). 
In order that i?^ preserves the symmetry T, we choose the reference metric 
g to also be a smooth Lorentzian metric which is stationary with respect to 
the vector field T, so that, 

g = -N{x) {dt + Ai{x)dx')'^ + hj{x)dx'dx^ (3.9) 

where we also assume here that T is globally timelike and bounded with 
respect to g so that iV > and bounded. Then hij gives a second Rie- 
mannian metric on the same manifold M.. Since now R^^ preserves the 
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stationary symmetry, the Ricci-DeTurck flow can be consistently truncated 
to a parabolic flow on the space of Lorentzian stationary metrics. Since this 
flow remains diffeomorphic to Ricci flow (subject at least to the normal com- 
ponent of ^ vanishing on any boundaries), we arrive at the interesting result 
that we may apply parabolic Ricci flow to stationary Lorentzian spacetimes. 
Likewise the Newton method will preserve the stationary symmetry and can 
be used to solve this elliptic stationary problem. 

In a situation where the solution we wish to find has a stationary Killing 
vector that is globally timelike and bounded then nearby to that solution 
the character of the Harmonic Einstein equations will be elliptic. Subject to 
imposing suitable boundary conditions on any boundaries or asymptotic re- 
gions, one may use the Lorentzian stationary Ricci-DeTurck flow or Newton 
method to solve for the solution. One must start with an initial guess that 
has globally timelike bounded stationary Killing field T, and then provided 
that guess is sufficiently good, one can hope the subsequent Ricci-DeTurck 
flow or Newton iterations preserve that T is globally timelike. 



3.3 Stationary black holes 

We now proceed to consider the case of non-extremal black holes. In the 
context of the discussion above now the norm of T will vanish either at the 
horizon itself, assuming that T is a globally timelike Killing vector (such 
as exist for certain Kerr-AdS black holes), or outside the horizon in the 
ergoregion. Since we are interested in the exterior of the horizon, in the first 
case we may treat the system described above for globally timelike T and 
now regard the horizon as a boundary where suitable boundary conditions 
are required. However in the latter, more general case, outside the horizon 
but inside the ergoregion we have the norm of T > and hence det hij < 0. 
Now the base manifold in the previous section fails to be Riemannian and 
then our argument above that the Harmonic Einstein equation is elliptic 
fails. 

In order to make progress we must use the Rigidity property of station- 
ary black holes, proved in D > 4 by Ishibashi, Hollands and Wald [35J for 
various asymptotics, including asymptotically fiat solutions. Assume there 
exists a stationary Killing vector T. Then the Rigidity theorem states that 
for an asymptotically fiat rotating black hole, so that T is not normal to the 
horizon, there exists a Killing vector K that commutes with T and which 
is normal to the horizon. Furthermore there exist some number A^ > 1 of 
commuting Killing vectors Ra, which also commute with T and asymptoti- 
cally generate spatial rotation with closed orbits of period 2tt. The theorem 
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states that K may be written in terms of these as, i^ = T + 0,"'Ra, for some 
constants il". Consequently the horizon rigidly moves with respect to the 
orbits of K in its exterior, and hence with respect to the asymptotic rota- 
tion generators Ra- Were this not the case one would expect gravitational 
radiation to be emitted from the region near the horizon and this would 
presumably violate the assumption of stationarity. 

Let us proceed by assuming Rigidity holds so that there exists a stationary 
Killing vector T and Killing vectors Ra for a = 1, . . . ,N, and the vector fields 
T and Ra are all commuting. We take the vectors Ra to generate spatial 
isometries with either compact or non-compact orbits. In the compact case 
we take the period to be 2it, and we allow axes of this symmetry where the 
isometry has fixed action. Rigidity implies we may write the normal K to 
our Killing horizon as, 

K = T + Q''Ra. (3.10) 

We may write the metric adapting coordinates to the isometries, 

ds^ = Gab{x) {dy"^ + Af{x)dx') {dy^ + Af {x)dx^) + hij{x)dx'dxi2,.ll) 

where y = {t, 2/"} and T = d/dt and Ra = d/dy"'. In analogy with the 
stationary case in the previous section we see that the geometry may be 
thought of as a fibration of the Killing vector directions over a base manifold 
M with metric hij . Technically Ai is the orbit space of the full Lorentzian 
spacetime with respect to the isometries T,Ra. We note that whilst in 4D 
for vacuum asymptotically flat solutions the circularity theorem implies one 
can find a coordinate system where the cross terms between base and fibre, 
Af-, vanish this is not expected to be the case for general stationary black 
holes in higher dimensions. At present, however, the only known solutions 
do in fact have vanishing Af. 

As with the analytic work on uniqueness, the aim now is to formulate the 
problem as an elliptic one on the orbit space M. We emphasize that here 
we are trying to find constructive numerical techniques to find black holes, 
rather than to prove their existence or uniqueness. With this in mind we 
make our key assumption; 

Assumption: (Ai, h) is a smooth Riemannian manifold. 

The full spacetime is Lorentzian, and so exterior to the horizon det (/^jy = 
det Gab det hij < 0. The chart breaks down at the horizon where det Gab = 
since the norm of K vanishes. It also breaks down at an axis of symmetry 
where some Ra vanishes and again det Gab = 0. However, our assumption 
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ensures that to the exterior of all horizons and axes then det Gab < and 
hence Gab is of Lorentzian signature. We regard the horizon and axes of 
symmetry of the RaS as boundaries for the base manifold M. We note our 
assumption above ensures that the geometry of these boundaries is smooth. 
For simplicity we assume here that the boundaries are only due to the hori- 
zon and vanishing of various Ra's. However more generally one might con- 
sider multiple Killing horizons, and boundaries where linear combinations 
of the RaS vanish J^^ I 

Harmark has discussed the above form of metric in the context of clas- 
sifying stationary spacetimes [3H]- The structure of A4 together with the 
data H."" at the horizon (or more generally horizons), and the data of which 
combination of RaS vanishes at the axis boundaries defines a 'rod struc- 
ture' for stationary spacetimes and has been conjectured to classify higher 
dimensional black holes. 

It is instructive to consider the simple example of the Kerr solution from 
this perspective of the time and rotation Killing directions being fibred over 
a smooth base. In the conventional Boyer-Lindquist coordinates the Kerr 
metric takes the form, 

ds'^ = Gttdt^ + 2Gt4,dtd(t) + G^^dcp"^ + hrrdr"^ + heodO^ (3.12) 



ith vanishing Af where, 




(A-a2sin2^) , ((r2 
Gtt = -^ ^ ^ , G^^ = sin2 0^^ 


+ a2)2 - Aa2 sin2 6) 




hoe = T, 1 


Ut(j) Ooin (7 , rirr . , 


ith A = r^ + 0^ - 2Mr and T, = r'^ + a^ cos^ 6. 


Here T = ^ and R 



(3.13) 

_ Jt 
~ d<j}- 

The outer horizon and axis are the boundaries of the base manifold M 
and are located at r = r^ (where A = 0) and 6 = 0, tt respectively. The 
Killing field K = T + QR is tangent to the horizon and timelike near there, 
where the angular velocity of the horizon is given as, O = 2? 2 • One finds, 
det Gab = —A, which vanishes at the horizon and axis, but not in their 
exterior. Whilst the 9 coordinate is a regular coordinate on the base at the 
rotation axes, the radial r coordinate is not at the horizon since A vanishes 
and so hrr — >• 00 there. We therefore define a new radial coordinate, p, such 
that dp = dr/vA and p = at the horizon, giving r = M+\/M'^ — a^ coshp, 
so that the components of the base metric hij are smooth at the horizon 
boundary. In particular in these coordinates the determinant of the base 

^^ This is discussed in detail in the case of D — 2 commuting Killing vectors I10| . 
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metric, 

hijdx'dx^ = -dr^ + ^dO^ = S {dp^ + dO'^) =^ det hij = T? >(3^14) 

and thus we see that since r/j > the base is indeed a smooth Rieman- 
nian manifold everywhere on and in the exterior of the horizon and axis of 
symmetry. 

3.4 Ellipticity of the stationary problem 

We note that we have not required the stationary Kihing field T to be time- 
like. In the presence of horizons it will become null on the horizon or be 
spacelike if the horizon is surrounded by an ergoregion. We reiterate that in 



the previous section ^3.2 it was precisely where T failed to be timelike that 
ellipticity would break down, since the base metric would fail to be Rieman- 
nian. The crucial observation is that for our class of stationary spacetimes 



(3.11), 



Rab = -Ig^^dad^gAB + ■■■ = -^/i™"5„9„GAB + . . . 
Rl = -Ig'^^dadpgAi + ... = -\h'^''drndn {GABAf) + ... (3.15) 

R?j = -\9''^Q<.dm + • • • = Af^'^dmdn [K, + GABAfA"^) + . . . 

where again the . . . represent lower than second order derivative terms. We 
see the equations have character determined solely by the metric hij , and by 
our assumption above that the base M is Riemannian, this is indeed elliptic. 
Ergo-regions may occur where T is no longer timelike, but our assumption 
that the base is Riemannian implies that some linear combination of the 
Killing directions T, Ka. is always timelike outside the horizons. 



In analogy with the previous section 3.2 , in order to ensure that -R^ shares 



the symmetries of g we choose the reference metric g so that T, Ra are again 
Killing with respect to it, and obey precisely the same assumptions as above 
for g. Thus we may write, 

ds^ = gi.vdX^dX'' (3.16) 

= Gab{x) {dy^ + Af{x)dx') {dy^ + Af{x)dx^) + hij{x)dx'dx^ 

and we further assume that {M, h) is a smooth Riemannian manifold. Then 
the Ricci-DeTurck fiow and Newton's method consistently truncate to the 



Lorentzian stationary spacetimes of the form (3.11) 



We must impose suitable boundary conditions at the boundaries of Ai. 
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Asymptotically we might impose the asymptotic flatness or Kaluza-Klein 
conditions mentioned above, which are compatible with ^ — t- 0. The new 
feature is that we have additional boundaries on the base corresponding to 
the Killing horizons and axes of symmetry and we will discuss this shortly. 
Using the Ricci-DeTurck flow or Newton method if we start from initial data 
in our stationary class, then for small flow times or updates we expect to 
remain in this class. In particular we expect {Ai, h) to remain a Riemannian 
manifold. Provided this condition holds for the solution of interest, and our 
initial guess is sufficiently close to this, then we might hope to reach this 
solution. 

We note that the maximum principle discussed in the static context relied 



on the inequality in (2.17), which results from the positivity of (V^^[/)(V^^'^) > 
for a Riemannian manifold. Since the equations for the static case are in- 
dependent of signature the maximum principle must equally apply in the 
static Lorentzian context. However, one can check that in the stationary case 
this term has indeflnite sign, and it is unclear if a maximum principle can be 
found ruling out solitons. This is not a problem in practice where one must 
simply check whether a solution obtained is a soliton or not. However, the 
elegant property that in certain static cases there can be no solitons does 
not obviously generalise to the stationary case. 



3.5 Boundary conditions for the stationary problem 

We conclude this chapter by now explicitly giving the boundary conditions 



for the metric components of our stationary spacetime (3.11) at the Killing 
horizon or symmetry axes. These are a generalisation of the boundary con- 
ditions determined for the 4D stationary axisymmetric vacuum problem of 
the classic uniqueness theorems [2]. They are derived and discussed in more 
detail in [31] and are consistent with the boundary conditions discussed by 
Harmark using particular coordinates on the base manifold [36j. Consider 
a Killing horizon with K = T + Q.'^Ra- It is then convenient to change 
coordinates as, 

t, y'^ -^ i=t, y'' = y''-nH (3.17) 

so that K = d/dt and Ra = d/dy"". Note that if y" is periodic then y"" is also 
a periodic coordinate with period 2-k. Now consider a boundary, either due 
to the vanishing of iT or a compact Ra- We take base coordinates x* = (r, x*) 
adapted to the boundary so that it lies at r = 0, and decompose the base 



36 Numerical construction of static and stationary black holes 

metric as, 

hijdx'dx^ = Ndr'^ + r N-.drdx' + h-.dx'dx^ (3.18) 

Horizon: For a Killing horizon we write the following metric components 
as, 

G,-^ = -r2/A, A^ = rg'', (3.19) 

for A = {i,y'') and then let X = |/a , 5^ , Gyayb , Af, N, Nj, h--\ be 
the set of functions describing our metric. Then by considering a change of 
coordinates, 

a = rcoshKt, 6 = rsinh«;t (3.20) 

and requiring in these Cartesian coordinates the components are smooth, 
we deduce the following behaviour is required in the stationary adapted 
chart; the functions X must be smooth functions of r^ and x* at r = 0, and 
furthermore obey a regularity condition, 

{f^-K^N)\r=o = (3.21) 

where k is constant and gives the surface gravity with respect to K. 

Axis: Consider an axis associated to a vanishing compact Ra- Without loss 
of generality choose this to be Rn- Then we choose to write, 

GyNA = r^fA , A^ = rg^ , (3.22) 

and let y = {/a , ff^ , Grt, %a , G-^-s , Af , N , Nj , h--} be the set of 
functions describing our metric (where a = 1, . . . ,N — 1). A similar analysis 
as for the shrinking of the Euclidean time circle previously implies that 
for a smooth metric we must have that the metric functions Y are smooth 
functions of r^ and x* at r = 0, and in addition we require, 

{fyN - N) \r=0 = (3.23) 

Of course we obtain analogous conditions for an axis with respect to a dif- 
ferent Ra- 
it is straightforward to check that the boundary conditions at the meeting 
of a horizon with an axis, or two axes, are compatible with each other. Take 
coordinates in the base x* = (ri,r2,x*) where ri = gives the position of 
the first boundary, and r2 = gives the second boundary, and hence the 
origin ri = r2 = is the meeting point. The boundary conditions near this 
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origin are simply the union of the boundary conditions for each boundary. 
Note this imphes that two boundaries (a horizon and axis, or two axes) meet 
in the base at right-angles. We reiterate that we have only considered axes 
arising from fixed points of the RaS, and more generally one could consider 
linear combinations of these vanishing. 

A very important point is that as discussed in the earlier section[2]l]having 
introduced boundary conditions we must check that these are compatible 
with finding Ricci flat solutions. To investigate this we must consider our 



choice of reference metric (3.16), which also is required to be regular and 
hence is subject to the same boundary conditions above for its components 
on the various horizon and axis boundaries. In particular we note that the 
surface gravity of the reference metric horizon must be the same as that of 
the actual metric. One can then explicitly check that, 

r|r-=0 = 0, dre\r=0 = 0, 9.e^|.=o = (3.24) 

both at a horizon and axis of symmetry, which is indeed consistent with 



the linear elliptic problem 2? • x = discussed in section 2.2 being well 
posed and admitting the trivial solution. Note that since ^^ = the Ricci- 
DeTurck flow should be diffeomorphic to Ricci flow in the presence of such 
boundaries. Furthermore the Harmonic Einstein tensor will be regular at 
the horizon and axis boundaries. Thus in our adapted coordinates it will 
also obey the same regularity conditions as the metric above. In particular, 
Ricci-DeTurck flow and Newton's method will preserve regularity, and will 
have the elegant result that they will leave the surface gravity constant. 

This work is dedicated to my father. 
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